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Bate  of  Chemical  Beacton  and  Its  Effeot  on  the  Process 
of  Droplet ' Combustion 

: 

7.  B.  Kuznetsov 

The  theory  of  combustion  of  a  fuel  droplet  in  an  oxidizer  stream  for  a  purely 
diffusion  case  has  been  developed  a  long  time  ago  [l,2].  For  the  general  case  is 
known  only  a  series  of  qualitative  results.  For  example,  it  was  established  that 
part  of  the  fuel  eyaporizing_  in  the_droplet_  doM  jnot_  combust  and  that  at  certain 

Fi’ST  l./vf  OP  T.'Tf.E 

conditions  the  droplet  itself  cannot  burn  L 3—5 J • 

Spaulding  made  an  effort  to  explain,  what  these  conditions  are,  but  the  base¬ 
less  disregard  of  convective  members  in  the  heat  conductivity  equation  led  to  the 
point,  that  the  value  V/R0  (ratio  of  rate  of  flow  coming  against  the  droplet  to 
its  radius)  during  flame  separation  (it  shows,  that  V/B  -function  only  of  the 
temperature  of  the  oncoming  flow)  is  almost  a  100  times  lower  than  the  experi¬ 
mentally  measured  by  Agafonova  [4]  and  Fedoseyev  [6]. 

In  experiments  [3.4]  v/ere  investigated  only  several  partial  instances,  by 
which  it  is  impossible  to  make  a  complete  picture. 

The  purpose  of  this  investigation  is  to  find  a  simple  method,  with  the  aid 
of  which  it  would  be  possible  to  evaluate  the  amount  of  uncombusted  fuel  and  to 
explain  the  conditions,  at  which  the  droplet  cannot  burn. 

Agafonova  [4]  derived  an  equation,  describing  the  process  of  combustion  of  a 
fuel  droplet  in  an  oxidizer  stream,  and  formulated  boundary  conditions.  To 
obtain  this  equation,  it  is  necessary  to  assume  that  1 

a)  the  boundary  layer  on  the  droplet  is  equivalent  to  a  certain  given  film, 
with  identical  thickness  everywhere,  ;  -  - 


b)  the  process  is  stationary,  c)  reaction  is  bymelecular  and  its  rate  is' 

diffusions!/ 

described  by  the  Arrhenius  equation,  d)  thermodiffusion  and/heat  conductivity 
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-  a  be  disregarded*  e)  const,  C^,ecn»t, 

*(*  +  «&  +  *l)(*  +  «W(>  +  W«)ip  (jjj  mg, 


This  equation  has  the  form  of 
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where 
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Here  T  -  temperature*  op  -  specific  heat  at  constant  pressure*  n^-concentia tioo. 
of  oxidizer  in  onooming  flow*  v-oass  rate  of  subs  tacos  within  the  film;  R-current 
radius  (read  from  the  center  of  the  droplet),*^ -  thickness  of  given  film;  D-diffusioa 
coefficient!  E-activation  energy j  K0-  pre-exponent*  q-thermal  effect  of  reaction 
by  oxidizer;  L  -  amount  of  oxidizer  in  g*neoessary  for  the  combustion  of  1  g  of  fuel; 
1*  speoifie  heat  of  evaporation;^-  density*  The  index  0  pertains  to  the  surface  of 
the  droplet*  1-  to  oncoming  flow* 

It  is  assumed*  that  f0-  boiling  point* 


Boundary  conditions  have  the  form  of 

•/ 


■i  * 


•it-v."*®** 

J9)  t 
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•  -r,;;£hejnagnitu4e'-X  in  the  interval  tyt)  1  ]  varies  slightly  in  comparison  with  the 
product  of  three  remaining  co-factors,  standing  in  the  right  part  of  equation  (l)» 
and  we  will  therefore  consider  that  X  »  const  *. 

In  equation  (l)  and  in  boundary  conditions  (3)  is  included  an  unknown  parameter 
uor^,  v/hich  cannot  be  determined  if  (l)  is  not  first  integrated.  In  reports  l4»5j 

p 

is  shown,  that  the  magnitude  u0r£  is  proportional  to  the  rate  of  evaporation, 

p 

practically  not  differing  from  (u0r‘)d.  (The  index  d  pertains  to  a  purely  iiffus- 

,ional  case)  In  conformity  with  this  we  will__assuae_,  jthat  ucr^  *>  (u0r^)d.  It  gives 

r;;,3-  l  ■;*  c « ?  r ::  _ ^  _±-  (4 ) 

*  1  — 

Prom  (3)  and  (4)  we  will  obtain 

I 

t5) 


We  will  introduce  a  value//}  =— * 


sL&\ 


dfiyst 


equalling  to  the  ratio  of  uncombusted 


L/rt  a /| 

fuel  flow  to  the  flow  of  evaporised  fuel. 


To  evaluate  same,  we  made  a  series  of  simplifications.  Their  validity  is 
discussed  below. 

Prom  (ij  and  (5)  we  have  » 

*U“— <«> 

We  shall  try  finding  the  value  of  the  integral  in  formula  (5). 

Since  in  the  reaction  zone  the  concentration  of  combustible  components  changes 
very  sharply,  and  the  temperature  to  a  lesser  ezter.t,  then  the  maximum  of  repro¬ 
ducing  the  concentration  and  maximum  rate  of  reaction  are  close  to  each  other. 

V.’e  wi:  1  a3cume  that  they  do  coincide.  It  gives 


c-o 


At  numorical  calculations  values  X  were  taker,  in  point 
pertains  to  maximum  rate  of  chemical  reaction.  ** 


The  index  a 
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Sj'vjr+b^  and  bj  with  an  accuracy  to  oonatant  multiples  art  equal 

to  oxidizer  tad  fuel  oonoantratlona  respectively. 

From  (7)  we  hair* 

*m -**»'*■*<  ,  (-) 

vbert 

2a, at  SthdLi&k  (5) 

•  *  ™"  I  •  ®  *“ 

«»  +  «t  «»+«». 

We  approximate  the  rate  of  reaction  (subintegral  function  in  (6})  by  &  jrad-al 
curve  to  at  it  was  made  by  Zel'dovich  and  Frenk-Kamenetskly  whan  solving  the  problem 
of  normal  flams  propagat  ionJjJ*  It  gives 

(0  [npx  v50 

*(♦)-'  g(*m.  «»)  »P*  Vi  <  $  <  V*. 

P  npa  ifc  <  <  1. 


Assuming  *  temperature  distribution*  corresponding  to  the  above  selected 

dependence  *g  (njr)  and  satisfying  the  boundary  conditions*  To  determine  tbeAljr  •  t|t 
vs  employ  an  idea  *  mentioned  in  the  above  given  report*  In  this  ease  it  leads 


to  s  relationship 


W  •  teemteuiwmy  01 

e  (^  — e*.  -  ®  —  — f  -  (5=; 

The  funot ion 'X is  approximated  by  two  straight  lines*  intersecting  ir.  point 

(  end  satisfying  the  boundary  conditions  (one  straight  lina  has  a  dip  -  a^» 

I  % •  ~^>n  ). 


and  the  other  one 


interral^ijr 


end  derivative 


♦ 

is 


easily  date  mined  by  these  straight 


lines t 


» W  «i-v/ 

a :  e. -'31 


|  ^  V—;  **  V*» 

Sines  ^  Tf  d-^Vl^  than  using  (8)  and  (10),  we  obtain  an  equation  for  the  deter* 

t  ~  . .  r  i  \ 

minstion  of  K  (vm -  &-«*)»  {:  -■ v«)  «*?  rr~ri  .  . 


*  (v«  -  W*  to*  -  4}»  a»  v-»  -  «*)*  U  -  vv.)  «?  rrrrj 
i  (flj  T  CjJ  (»Vm  V  “ — 
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Hiring  found  frcm  (12),  6^  from  (8)  ind  (12),  we  will  obtain from  (11)* 
Next  in  found  cj^. 

Suob  calculation*  wn  mad#  for  an  alcohol  droplet  in  an  air  a  tram  (table  )• 

Zn  flg.l  i a  given  the  relationship /^(E)  at  const.  We  investigate  behavior  of 
functions'*^  (K)  near  the  point  Since  the  raluea  a  +d«  1*  o  'jr**  ^ 

do  not  change  muoh»  for  an  approximate  analysis  they  can  be  considered  constant*  Equa¬ 
tion  (12)  acquires  the  form  of 


'  const. 


/12a) 


At  K-^  and  b  -  const  (it  corresponds  to  a  rise  in  rat*  of  chemical  reaction 

to  infinity  or  a-  reduotion  to  zero  in  the  rate,  of  flow*directed  around  the  droplet* 

or  stoichiometric  coefficient)^ — ^  tlie  ileat  generating  ability  of  the  fuel 

or  concentration  of  oxidizer  in  the  stream  increases*  K  decrease# «but  1/b  exp  (b/ctr 

(IT] 

-fd)  rises  much  sharper  and-vj^  also  tends  toward 


?ig*2* 

Since  ^  <*> 

and  the  amount  of  uncocbusted  fuel  tends 


to  zero. 
At 


f %2  ^4.  bp  ^0  (negative  concentration  of  fuel),  and  at  -  d/e* 

0#  consequently  -  d/04 ^  4  \|fcd*^y 

Ve  will  now  examine  the  conditions  /  flame  separation.  Ve  will  mention  first  of  all 
that  at  greater  values  of  the  Nuseelt  criterion  R0/V  Jit  is  asstned*  that  HU  - 
2-j-0*6*V^®")*  Consequently  at^nchanged temperature  of  oncating  flow  flame  separation 
will  take  place  at  one  and  the  very  same  3^/V  values. 
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It  ve  increase  the  rat*  of  flov  V.tben  tb*  rat*  of  evaporation  will  be  increased. 
In  the  and.  tbere  abould  b#  flan*  separation, 
la  the  interval  d/«(  the  ▼»!*»* 


U 


(12'=) 


1*  limited  from  the  top*  consequently  at  oertain  K  values  equation  (12)  baa  no  aolu  — 
tione .  The  latter  mean*  that  stable  combustion  of  the  droplet  is  impossible*  The 
critioal  value  K  is  found  from  formula 


dK 

** m 


<0. 


f'l) 


which  determines  a  certain  curve  which  appears  to  be  the  boundary  between 

the  zone  • where  stable  combustion  is  possible •  and  the  zone*  where  stable  combustion 
is  impossible*  Such  a  curve*  calculated  for  a  drop  of  alcohol  ^burning  in  the  flow 

of  air*  is  presented  in  fig*2.  . 

of/ 

To  cheek  the  accuracy/  the  approximate  method  of  solving  equation  (1)  was  made 
a  comparison  of  val^  d 6e  /d>|r^  found  from  an  accurate  solution  by  numerical  in 
tegration  and 


by  the  above  explained  method.  Differences  have  not  exeeeded  10- 


13%. 
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Fig*3  Fig.4 

Calculated  were  also  the  dependenoea  of  the  amount  of  un combusted  fuel  (gasoline) 
upon  the  rate  of  oncoming  air  stream  and  the  relationships  between  the  temperature 
of  the  flow  and  its  rata  during  flama  separation. 

The  results  of  calculations  are  shown  in  fig«3  and  4  (dotted  curves).  For  com¬ 
parison  on  the  very  a  am*  figures  are  given  experimental  curves  borrowed  from  r4~  • 
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Numerical  values  of  magnitudes,  necessary  for  the  calculation,  are  given  in  table 

Table 


Name 

E.cal/mol 

L 

1 

i  i 

q.  A* a/4  4,  KCJ/i 

1  1 

yt/Mp 

t 

i  .  I 
i  *  i 

r,\K  | 

I  1 

I 

!  *. 

D  j)  ,  g. cm/sec 

Ethyl  alco¬ 
hol 

Gasoline 

1 

1 

\  42500 

1  45000 

2,1 

3,43 

3*00 

3000 

220 

70 

0,33 

0,33 

i 

0.23 

0.23 

340 
365  | 

7-10»» 

7-10“ 

3.2-10-* 

1,72-10 r* 

1 

Greater  discrepancies  between  values  obtained  experimentally  and  theoretical 
ly,  are  explained  by  failure  to  consider  the  afterburning  along  the  path  after  the 
droplet. 

The  calculated  dependence  of  the  rate  of  the  flow  during  flame  separation 
upon  temperature  of  the  flow  is  in  excellent  conformity  with  experimental  data. 

Between  the  values  V/Hq  during  flame  separation,  obtained  experimentally  ar.d 
theoretically,  there  is  also  satisfactory  conformity  (at  370°  they  constitute  pSOO 
and  5^00  sev  ), 

The  author  expresses  sincere  thanks  to  L.  A.  Xlyaohko  and  M.  P.  Samozvantsev 
for  critical  remarks  and  valuable  council. 
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